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A Singularity Problem with f(R) Dark Energy 
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In this paper, I point out that there is a curvature singularity problem appearing on non-linear 
level that generally plagues f(R) models that modify Einstein gravity in the infrared. It is caused 
by the fact that for the effective scalar degree of freedom, the curvature singularity is at a finite 
field value and energy level, and can be easily accessed by the field dynamics in the presence of 
matter. This problem is invisible in linearized analysis, except for the tell-tale growing oscillatory 
modes it causes. In view of this, viability of many f(R) models in current literature will have to be 
re-evaluated. 
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What is causing the observed accelerated expansion of 
the Universe today is one of the biggest open question 
in modern cosmology. Trying to explain it by modifying 
theory of gravity rather than by introducing a mysteri- 
ous dark energy component has been a popular pursuit 
as of late. Unfortunately, it is proving to be a rather dif- 
ficult thing to do consistently, while avoiding variety of 
stringent observational tests of gravity we have at our dis- 
posal. A class of such models that received much atten- 
tion recently is the one which modifies Einstcin-Hilbcrt 
gravitational action by replacing Ricci curvature scalar 
by an arbitrary function of the curvature 
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Introduced in cosmological context for the case which 
modifies gravity in the high energy limit in a seminal pa- 
per by Starobinsky [l[ and studied in [1, 0, Hj] , this model 
has later been adopted for infrared modifications of grav- 
ity as well [f| H[ . For the latter application, it turned out 
to be not without problems. Certain constraints have to 
be imposed on function f(R) for the model to be linearly 
stable Q and cosmologically viable d, H, [l(| • The first 
attempts failed these constraints right away, but since 
then, models that evade them have been found (for ex- 
ample see [III [TH, [l3| and references therein) and enough 
trust has been placed in their viability to study cosmo- 
logical structure formation in detail [lj, EH • 

In this paper, I point out a serious curvature singu- 
larity problem that affects many, if not all, infrared- 
modified f(R) models. Being non-linear in nature, it 
escaped scrutiny so far. 

As it is well known, a new scalar degree of freedom ap- 
pears in f(R) gravity that is not there in Einstein theory 
(sometimes dubbed the scalaron). Conformal transfor- 
mation of the metric can be employed to make it explicit 
in the action [1, 0|- In this paper, I will avoid doing that 
to keep the usual matter coupling to the metric, and 



work with the action |T]) directly. Variation with respect 
to metric yields gravitational equations of motion 

f'Rpu - f! tfl „ + ( □/' -~fj 9llv = SttG t m „, (2) 

where prime (') denotes the derivative of the function 
/ with respect to its argument R, and □ is the usual 
notation for covariant DAlembert operator □ = V a V a . 
The equation of motion for a new scalar degree of freedom 
is given by the trace of equation ([2]) 



□/' = I (2/ 
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Identifying the scalar degree of freedom explicitly by a 
variable redefinition 



(4) 



the equation ([3]) above is cast in the form of equation of 
motion of a canonical dimensionless scalar field 4> with a 
potential V and a force term T 

U<$> = V\$)-T. (5) 

The effective scalar field potential V{<f>) is determined by 

V'(<f>) 



(6) 



expressed in terms of the scalar variable <f>. In practice, 
given f(R), it is usually difficult to invert the definition 
of the scalar degree of freedom ((4]) explicitly, so it might 
be more convenient to determine effective potential V in 
a parametric form instead. By integrating 



dR~ dchdR- 3 (2/ f ' R)f " 



(7) 



potential V((f>) is then given by a pair of functions 
{<j)(R),V(R)}. The force term T that drives the scalar 
field (j) is a trace of the stress-energy tensor T, which for 
perfect fluid is simply 
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Let us consider a homogeneous cosmological model 
in f(R) gravity, with the usual complement of matter 
fields. Expansion of the Universe is described by a flat 
Friedman-Robertson- Walker metric 



ds 2 = -dt 2 



\t) dx 2 



(9) 



and the scalar gravitational degree of freedom <f> obeys a 
usual scalar field equation, albeit with a force term on 
the right hand side 



ZH(j> + V'{(j)) =T. 
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The analog of Friedman equation in f{R) cosmology is 
not so transparent. Let us consider tt component of grav- 
itational equations of motion @. For metric ([5]), it is 

3H(f'y-3-f' + y = 8wGp. (11) 
a 2 

Note that unlike the usual Friedman equation, higher 
derivatives of scale factor a appear. Second derivative d 
is written out explicitly, and a third derivative is hiding 
in a time derivative of /' term, which itself contains Ricci 
curvature, and hence d. Seeing a second derivative of the 
scale factor, one might be tempted to treat the above 
equation (llip as a dynamical evolution equation for the 
scale factor. Doing so, however, is not a very good idea. 
For small deviations from Einstein gravity, the coefficient 
in front of d goes degenerate, and the equation (fTTjl does 
not have a good limit determining d (which is not all 
that surprising, considering that Friedman equation in 
Einstein gravity does not constrain d directly) . To get a 
proper limit, let us instead get rid of d in favour of the 
curvature scalar 



, a a 

R = 6 I - + 

a a 



After that is done, the equation (fTTj) becomes 



H 2 + (\nf'YH 
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and its role as a constraint equation is revealed. In the 
limit of Einstein gravity /' — ► 1, and so the last two terms 
on the left hand side disappear, and one is left with the 
usual Friedman equation. In the general case, the extra 
terms are functions of scalar degree of freedom (f> and its 
first time derivative. No higher derivatives appear in this 
equation anymore. 

Thus the following simple picture of dynamics in the 
f(R) cosmology emerges. Above the infrared modifica- 
tion scale i?0) the expansion rate of the Universe is set 
primarily by the matter density, just like in the usual cos- 
mology, with small corrections. Only once the local cur- 
vature drops below Rq, the expansion rate starts feeling 
the effect of gravity modification. The spacetime curva- 
ture, on the other hand, is controlled by the scalar degree 
of freedom <f> which gravity acquires. It obeys the usual 
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FIG. 1: Effective potential of a scalar degree of freedom in 
f(R) gravity model (|14p with A = 2 and n = 1. Diamonds 
mark the location of critical points. The part relevant to 
cosmological evolution is emphasized by thick blue line. 



scalar field equation (flTJf with potential V(<p), the shape 
of which is directly determined by function f(R), and a 
driving term from the trace of matter stress-energy ten- 
sor. 

But here is the problem: it turns out that precisely 
those functions f(R) that lead to Einstein-like gravity 
action in the large curvature regime, yield a potential V 
with an unprotected curvature singularity. 

As a case in point, consider Starobinsky's disappearing 
cosmological constant model [l2 |. which has been very 
carefully constructed, and avoids all known linear insta- 
bilities. It is described by 



f(R) = R + X (l + R 
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(14) 



where I have taken a liberty to absorb the cross-over 
curvature scale Rq into rescaling of coordinates (which 
become dimensionless and are measured in length units 
corresponding to Ro)- For definiteness, let us take n = 1. 
The scalar degree of freedom in this model is given by 
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in terms of curvature, so large curvature limit R — > ±oo 
corresponds to <j> — 0. Flat spacetime with R — also 
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FIG. 2: Adding matter destabilizes the vacuum. Although 
effective potential inside constant density matter distribution 
still has a minimum, it is very shallow, and cannot protect the 
field <j> from reaching curvature singularity X, which becomes 
energetically accessible from asymptotic vacuum state B. 



corresponds to <$ = 0, which gives us a hint that the 
potential is going to be a multi-valued function. The 
potential can be evaluated by integrating ([7]); up to an 
arbitrary constant it is 



V = 



\ 2 R{3 + IIR 2 + 21i? 4 - 3i? 6 ) 

24(1 + i? 2 ) 4 
Ai? 2 (l + R 2 - R A - R 6 ) 



(16) 
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The effective scalar potential is plotted in Figure [T] for 
A = 2, and is indeed multivalued. Let us walk through 
the interesting locations on this plot. Point A is a posi- 
tive curvature singularity R = +oo. Point B is the stable 
de Sitter minimum in this model, and point C is the un- 
stable de Sitter maximum; their curvatures depend on A. 
Point E corresponds to a flat spacetime, which although 
a solution in this model, is unstable. Points D and F are 
critical points with /" = that occur at R = ±1/V3; 
potential branches there. Finally, point G is a negative 
curvature singularity R = — oo. Only the small part of 
this potential is actually relevant for cosmological evolu- 
tion from initial singularity to today, and it lies in the 
arc AB, shaded blue in the Figured] 

The most striking feature of the potential in Figure [1] 
and the core of the problem for infrared-modified f(R) 
models, is that curvature singularity at point A is finite 
distance away both in field and energy values from the 
place we are supposed to live in. Scalar degree of free- 
dom (j) directly feels the matter distribution through the 
force term; for equation of state w < 1/3 the force is di- 



rected to the right, and drives the field <fc up the wall to- 
ward point A and infinite curvature. Characteristic scale 
of the potential V is the cross-over curvature scale Ro, 
and hence of the same order of magnitude as a present 
day cosmological constant, which is exceedingly low com- 
pared to matter densities we encounter every day. Given 
the scales involved, it appears to be quite easy to over- 
drive the scalar degree of freedom and make it "jump 
out" of the potential well by doing simple manipulations 
with normal matter (say a pile of dust), which would 
cause catastrophic curvature singularity. Needless to say, 
if this were to happen, it would not make for a desirable 
(or even viable) model. Similarly, but less dramatically, 
matter with sufficiently stiff equation of state can desta- 
bilize the model by driving the field to the left past the 
unstable point C . 

The presence of the curvature singularity a finite dis- 
tance away is extremely disturbing by itself, but let us 
examine more carefully if it is reached by physically rea- 
sonable solutions. Inside a constant density matter dis- 
tribution, one can think of a (constant) force term T as 
coming from a linear field potential T{(\>* — 4>) instead, 
and introduce a new "in matter" effective potential 



U(<t>)=V(<l>)+F(cj>* 



(17) 



where denotes the asymptotic de Sitter vacuum field. 
The comparison between the two potentials V and U for 
T = Rq is shown in the Figure[2] As you can see, addition 
of matter slopes the potential U, shifts the (stable) min- 
imum to the right but makes it more shallow, and lowers 
the curvature singularity point X. The density needed 
to make the curvature singularity energetically accessible 
from vacuum B is given by the ratio of potential barrier 
5V = Va — Vb to the scalar field value distance from 
vacuum to singularity 5<j) = 4>a — 4>b = — </>*■ It is of the 
same order as the density of dark energy today, with a 
numerical factor which depends on the model. So for vast 
majority of physical solutions with matter, the curvature 
singularity is energetically accessible from asymptotic de 
Sitter vacuum, and the potential minimum is so close to 
curvature singularity that it would be invisible if actually 
plotted to scale. 

Energetical accessibility of the curvature singularity 
causes problems. For example, if one takes a cosmologi- 
cal solution approaching dark energy domination today, 
and traces it back into the past, one is very likely to en- 
counter a curvature singularity. This has been noticed 
numerically [H, but the underlying reasons for it 
and the extent of the damage were not fully realized. It 
is also most likely the cause of growing oscillatory cur- 
vature modes [12j which signal the break-down of linear 
expansion due to closeness of potential minimum to cur- 
vature singularity. Although a more detailed analysis of 
the approach to singularity is in order, from equation 
(fl~3"|) it appears that the singularity occurs at finite red- 
shift, density and expansion rate, and is driven by the 
divergence of the second derivative of the scale factor 
a (which would make it rather weak, but a singularity 
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nonetheless). 

Although I focused on cosmology so far, perhaps a 
more deadly argument against having a curvature sin- 
gularity at finite distance in the field space comes from 
considering a gravitational field of a static dense com- 
pact object (like a neutron star). Although the exact 
non-linear solution of this problem is more complicated 
to analyze [H, [l9| and is beyond the scope of this article, 
I can give a very simple estimate if the problem occurs. 
As we have seen, the energetics of the scalar gravitational 
degree of freedom are by far dominated by the matter 
driving term. If we discard the contribution of non-linear 
potential V (which is negligible everywhere except maybe 
very close to singularity at <f> = for compact object), the 
equation for gravitational field of a static matter distri- 
butions becomes a simple Laplace equation 



A0 
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(18) 



Comparing this with an equation for Newtonian gravita- 
tional potential 



A$ = AirGp, 



(19) 



we get a simple estimate for the excitation of scalar grav- 
itational degree of freedom in f(R) gravity in terms of 
Newtonian potential well depth $ of the compact object 



(20) 



where 4>„ is the asymptotic value of a t infinity, i.e. 
the minimum value 4>b — —S<f>. But unlike Newtonian 
potential $, which has to diverge to cause singularity, 
or reach —1/2 to form a horizon, gravitational degree of 
freedom cj) needs only to change by a (small) amount 8<j) 
from its vacuum value to create a singularity. So unless 
an infrared-modified f(R) model leads to a potential with 
curvature singularity separated from vacuum by at least 
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one would end up with a curvature singularity without 
horizon in a compact astrophysical object like a neutron 
star. This condition is rather easy to violate unless spe- 
cial care is taken in model-building. For example, for 
Starobinsky's model (fl4"|) with n = 1 and A = 2 (as in 
Figure [T]) 8<fi — 0.0874 <C 1/3, and is even smaller for 
larger values of A, for which it decreases as 5(j) ~ (2A) -2 . 
Since in general one needs /' > for graviton not to 
be a ghost, one would need —1 < 0* < —1/3 to avoid 
both problems, the prospects of achieving which without 
fine-tuning do not look good. 

This curvature singularity problem is in no way 
unique to Starobinsky's disappearing cosmological con- 
stant model 12j , which I have taken as an example simply 
because it is one of the most carefully constructed models 
so far. In fact, any infrared-modified f(R) gravity model 
suffers from it. Let us consider arbitrary function f(R), 



and require that it reduces to Einstein gravity for large 
curvature, and has an analytic expansion 



f(R) = r + a + -L y ^ 

v ' R a ^ R n 



(22) 



with a leading term no/R a (with a > 0). Then the 
leading terms for large R asymptotic behavior of scalar 
gravitational degree of freedom (J4| and potential dTj) are 

, «M0 -rr , (tt + l)j«0 /f)Q x 

<b ~ -r, V — const . (23) 

Y R a + l 3R a v ' 

The value of <fi goes to zero in large curvature limit, and 
the potential V has power law dependence on <j> 



v(4>) 



const — 



(a + l)flg 
3 \ap \ 



7 



101" 



with exponent 7 valued between zero and one 

a 



7 



a + 1 



(24) 



(25) 



Thus, the values of both the field and the potential at cur- 
vature singularity are finite for a generic f(R) infrared 
modification of gravity which recovers Einstein gravity 
perturbatively in the large curvature limit. This means 
the arguments I made above apply generically, and via- 
bility of many f(R) models in current literature will have 
to be re-evaluated. At the very least, the bound (f2"Tj) will 
have to be satisfied for the model not to be ruled out 
immediately. But even if the estimate for compact ob- 
jects I made here looks OK, any infrared-modified f(R) 
models should be scrutinized very closely for dangerous 
curvature singularities that could be present. In a sense, 
infrared modification of f(R) gravity forces one to con- 
front the question of ultraviolet completion of the theory. 

Finally, let me comment on how this problem looks 
like in equivalent scalar-tensor theory formulation 0, 0] • 
Conformal transformation to an Einstein frame with met- 
ric ds 2 = /' ds 2 turns the scalar degree of freedom into a 
canonically normalized scalar field ip with potential 



t 



fin/', 



1 



W{^) = -e-^{Rf - f). (26) 



The asymptotics of scalar degree of freedom in Einstein 
frame are very similar to the above story: the field ip goes 
to zero in large curvature limit, and the potential has the 
same unprotected power law asymptotic W ~ a — b\ip\' y . 
But where did the singularity go? The answer is subtle: 
while the conformal factor itself appears to be regular 
(/' — » 1), its second derivatives are not (potential deriva- 
tive blows up as I tp | 7-1 in equation of motion), which 
can cause a curvature singularity in Jordan frame even if 
Einstein frame metric was regular. 
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